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ABSTRACT 


Cofty 


The  equations  of  gas  dynamics  in  a  tube  with  varying  cross  section  are  an 
example  of  a  nonhomogeneous  system  of  conservation  laws.  In  this  work  we 
study  the  Riemann  problem  for  this  system  by  viewing  it  as  a  perturbation  of 
the  classical  equations  of  gas  dynamics  in  a  uniform  tube.  We  also  study  the 
Riemann  problem  and  the  formation  of  singularities  for  a  related,  but  simpler, 
problem  of  a  nonhomogeneous  Berger's  equation. 
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significance:  and  explanation 


A  one-dimensional  mathematical  model  for  the  flow  of  qas  in  a  nozzle  (a 
tube  with  varyinq  cross  section)  takes  the  form  of  a  nonhomoqeneoiis  system  of 
conservation  laws.  It  is  an  interestinq  problem  to  study  how  closely  this 
one-dimensional  system  models  the  real  two-  and  three-dimensional  problem. 
With  this  qoal  in  mind  we  study  two  problems:  a)  the  Riemann  problem  (i.e., 
the  solution  to  the  problem  whose  initial  datum  is  a  step  function)  for  the 
one-dimensional  system/  b)  the  formation  of  shocks  for  a  sinqle  first  order 
equation  which  has  most  of  the  structure  of  the  one-dimensional  model  system. 


The  responsilDility  for  the  wordinq  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC/  and  not  with  the  a\ithors  of  this  report. 
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Fhao-Shiunn  Lin  an<^  Rp?ra  Malek-Madani 

1*  mtrnHnctlon.  The  equations  of  qas  dynamics  in  a  uniform  tube  have  been  studied  auine 
extensively  in  rerent  vears.  It  is  well  known  that,  as  a  hv{>erholic  conservation  law,  these 
equations  exhibit  discontinuous  solutions,  while  the  initial  value  problem  is  not  mathemat¬ 
ically  well  posed  in  the  class  of  weak  solutions  f1|*  It  is  not  difficult  to  envisaqe  the 
mathematical  reason  for  the  nonsmoothness  of  solutions.  These  equations  enjoy  a  full  set 
of  real  characteristics  and,  if  the  initial  values  are  chosen  proper IVr  the  information 
carried  by  the  characteristics  will  overlap  and  shocks  develop.  The  problem  under  study  in 
this  paper  has  one  additional  property,  namely  the  variation  in  the  tube's  cross  section, 
that  will  persumahly  contribute  even  further  to  the  shock  producinq  mechanisms. 

Section  2  concerns  with  the  derivation  of  the  equations  studied  in  this  work.  The 
arguments  of  Hughes  f21  have  been  followed  and,  as  it  will  become  apparent,  the  system 
under  consideration  is  an  example  of  nonhomogeneous  hyperbolic  (Tonservation  laws.  In 
Section  3  a  simpler  but  related  problem  is  discussed  for  the  purpose  of  understandina  the 
shock  producing  mechanisms  that  do  not  exist  in  the  homogeneous  problem# 

Section  4  concerns  the  solution  of  the  Riemann  problem.  It  is  well  known  (31,  [41 
that  the  solution  of  the  Riemann  problem  played  an  essential  role  in  developing  a  numerical 
scheme  in  order  to  solve  the  initial  value  problem  for  the  equations  of  aas  dynamics  in  a 
tube  with  uniform  cross  section.  Motivated  by  this  fact  T.  P,  Lui  [51  applied  a  modified 

Riemann  problem  for  the  general  order  nonhomogeneous  conservation  laws  and  developed 

an  iterative  scheme  which  converges  to  the  weak  solution  of  the  initial  value  problem. 
Although  the  above  scheme  is  guite  successful  theoretically  it  is  rather  difficult  to 
implement  it.  Since  we  have  in  mind  a  concrete  example  from  the  equations  of  gas  dynamics 

it  is  our  contention  to  propose  a  simpler  Riemann  problem  and  hope  that  it  would  give  rise 

to  more  manageable  computations.  We  are  presently  studying  this  problem. 


Sponsored  hv  ♦'he  United  States  Arrv  tinder  f’optrart  No,  nAAG2^-Hn-\ ’ -r^'U i  ,  This 
ma^er^al  is  based  upon  work  sur'tv^rt-ed  bv  the  National  Fcienre  Fntir.ht  i- 1  nn 
Grant  No.  MrS-7U27062. 


2.  rn^rivation  of  th»»  modPl  pouation.  Consider  an  inviscid  isentropic  oas  flow  throiiqh  a 
two  dimensional  duct  V  -  (  (x,v)|a^(x)  <  y  <  A^fx),  r  x  <  "*} .  lyie  motion  of  the  qas  is 
ooverned  by  the  equations  of  conservation  of  mass  and  linear  momentum 

p  ^  +  (pu(  t  =  0  > 

(Pu).  +  (Pu^  +  P)  +  (puv)  =  P  ,  (2.1) 

t  XV 

(pv)^  +  (puv)  +  (Pv^  +  p)  =0 

t  X  V 

with  p  =  f(p  ),  where  P  =P(x.y,t)  is  the  density,  p  =  p(x,v,t)  is  the  pressure  and 
u  =  (u,v)  is  the  velocity  vector,  toqether  with  the  Neumann  boundary  conditions 
u(x,A^  (x)  ,t)AMx)  =  v(x,A^  (x),t),  i  =  1,2 

and  the  initial  conditions 

P  (x,y,0)  =  P|^(x,y)  , 
u(x,y,0)  =  UQ(x,y)  , 
v(x,y,0)  =  v^(x,y)  . 

In  the  remainder  of  this  section  we  will  outline  briefly  the  procedure  discussed  in 
(21  which  approximates  (2.1)  by  a  one-dimensional  nonhomoqeneous  system  in  the  variables 
p  and  u.  For  a  physical  quantity  q(x,y,t)  defined  in  the  reqion  V  we  define  the 
averaqe  (  q  >  of  q  in  the  y-dlrection 

,  A2(x) 

*  =  TTTT  /  a(x,v,t)dv 

'  ’  A^(x) 

where  A(x)  =  A2(x)  -  A^(x).  Averaqinq  each  equation  in  (2,1)  and  usinq  the  boundary 
conditions  yield 

A*  (x ) 

<P>t^<P">x=  <*’«>' 

(pu>^  t  <pu^^  *  (P>^  =  <pu^,  (2.2) 

A*  (x  ) 

(pv>^  +  (puv)^  +  (Py>  =  -  (PUV>  , 

(  P)  =  (f  (p))  . 

In  order  to  further  simplify  (2.2)  we  make  the  followinq  assumptions; 
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(A) 

the  total  variations 

of  A^(x) 

and  A2(x) 

arp  snail. 

(B) 

the  quantity  |^t  << 

1 1  i f 

th#»  flow  is 

t>r<*domir>ant ly  in  y-^irpction 

(C) 

(f  (P  )  >  =  f  ((p  > )  for 

sotne  f  • 

Then  it  is  reasonable  to  assune  that 

<  Pu>  =  <  P><u)  , 

2  2 
(  PU  )  =  <P)  (U) 

etc.  An  asymptotic  analysis  with  respect  to  l-jl  adds  more  plausihility  to  the  conations 
(2.3).  Thus  (2.2)  becomes 

A*  f  Y  ) 

%  ^  • 

(PU)^  +  (PU^  +  p)^  =  -  pu^  ,  (2.  3) 

p  =  f(p ) 

where  we  have  made  the  followinq  identifications 

<p)  ~P<>'.t).  (u>~«(x,t) 

etc. 

System  (2.3)  is  the  one-dimensional  approximation  of  (2.1).  It  should  be  pointed  out 
that  as  far  as  the  authors  )<now  there  has  not  been  a  riqorous  analysis  of  how  reasonable 
the  assumptions  (A-D)  are.  Nevertheless,  the  system  (2.3)  is  a  mathematically  tractible 
model  of  (2.1).  It  is  believed  that  the  study  of  (2.3)  will  shed  some  liqht  to  the 
structure  of  the  solutions  of  the  more  difficult,  but  exact,  equations  of  qas  flow  in  a 
duct  with  variable  cross  section. 
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3.  Formation  of  singularities  for  the  equation  =  a(x)u.  “efore  proc<’o3 l r.g 

with  the  solution  to  the  Riemann  problem  for  the  system  (2.3)  it  is  instructive  '  >  study 
how  the  spatial  dependence  of  (2.3)  enters  as  an  important  feature  in  producing  shoc)cs. 
The  nonhomogeneous  Burger's  equation 


u  *  ^(u)  =  a  (x)u 

t  X 


(3.1  ) 


u(x,0)  =  u  . 

has  most  of  the  essential  structure  of  equations  (2.3)  with  one  additional  advantage  that 
(3.1)  has  only  one  family  of  characteristics.  Our  goal  in  this  section  is  to  show  that  the 
global  smoothness  of  the  solutions  of  (3.1)  depend  strongly  on  the  sign  of  a(x)  and 
a'(x).  By  way  of  example,  consider  the  simple  case  of  =V2“^*  I*  a(x)  was  a 

constant  a,  then  (3.1)  can  be  solved  explicitly  along  characteristics  to  give  us  a  smooth 
solution  along  a  family  of  parallel  characteristics.  The  slope  of  these  characteristics 
depends  on  a  and,  as  easily  seen,  it  is  an  increasing  function  of  a.  Therefore,  it 
seems  plausible  that  as  a(x)  decreases  the  characteristics  will  intersect  in  finite  time. 

The  definition  of  a  shocic  discussed  in  the  theory  of  conservation  laws  lends  itself 
naturally  to  equation  (3.1).  Let  a  family  of  characteristics  x(t,5)  be  defined  by 

x(0,5)  =  i  . 

We  assume  throughout  this  paper  that  equation  (3.1)  is  hyperbolic  and  genuinely  nonlinear, 
i.e.,  ^'(u)  >  0  and  ^"(u)  >  0  respectively.  Under  these  conditions,  a  singularity 

(shock)  develops  if  two  characteristics  intersect.  Let  x(t,5^>  and  x(t,52)  t*  two 

characteristics  Intersecting  at  {T,x*),  that  is,  x*  =  x(T,5^ )  =  Further, 

dx(T,t^) 

suppose  that  the  two  characteristics  actually  cross  each  other  and  - r: -  * 


II  -  ^•(u(x,t)) 


dt  dt 

It  follows  from  the  above  definition  of  characteristic  curve  that 
'^'(u(x(T  ,5^),T  ))  *  ^'(u(x(T  )).  Since  is  assumed  to  be  monotone  the 

solution  u(x,t)  becomes  multivalued  at  (T,x*)  and  a  shoe)!  develops.  It  should  be  noted 
that  a  shock  formed  due  to  the  interaction  of  two  characteristics  in  forward  time  already 


has  the  entropy  inequality  embodied  in  it  (cf.  (1)) 


An  indication  that  the  smooth  solution  u(x,t)  will  be  discontinuous  at  a  point 
(x,t)  IS  that  Uj^(x,t)  becomes  unbounded  in  finite  time  [6].  On  the  other  hand,  u^^ 
can  be  evaluated  along  the  characteristic  x(t,C)  as 


u  (x!t,C),t) 

X 


U^(x(t,C),t) 

x^(t,?) 


Thus  u  will  become  unbounded  if  x,(t,C)  approaches  zero.  This  presents  an  alternative 

X  ^ 

way  of  establishing  the  formation  of  a  shock  (c£.  (71). 

The  following  lemma  states  that  if  x^(t,^)  becomes  zero  in  finite  time  then  two 
characteristics  must  intersect.  Let  v(t,C)  5  x^(t,5).  Note  that  v(0,C)  =  1. 

Lemma  3. I :  Suppose  that  there  exists  T  <  “  and  5  such  that 

v(T,4)  <  0  . 

Then  there  are  5^  and  with  x(T,C^)  =  Moreover,  if  u  satisfies  an 

equation  of  the  form  g(x/U)  along  the  characteristic,  it  follows  that  the  solution 

develops  a  shock  in  finite  time. 

Proof ;  By  way  of  contradiction,  suppose  that  for  all  5^  *  ^2' 
implies  that  the  function  f(C)  defined  by 

£(5)  =  x(T,?) 

is  monotone.  Therefore,  f'(C)  will  always  be  nonnegative  which  contradicts  the 
hypothesis.  Hence  there  are  two  characteristics  and  which  meet  at  (x,T)  for 

some  X.  On  the  other  hand,  by  the  standard  uniqueness  theorem  in  ordinary  differential 
equations,  the  above  characteristics  viewed  in  the  (x,u>  plane  reach  the  line  x  =  x  at 


two  different 

values  of 

u 

at  time  T. 

,  This 

completes  the  proof  of  the  Lemma. 

Theorem 

(3.1); 

;  Suppose 

a  C 

c’  10,-), 

Uq  is 

a  position  constant,  d'(u)  >  0 

k  >  0 

for 

some 

k: 

(>) 

If 

a(x) 

<  0 

then 

(3.1)  has 

global 

smooth  solutions. 

(2) 

If 

a(x) 

>  0 

and 

a'(x)  >  0 

then 

(3.1)  has  global  smooth  solutions. 

(3) 

If 

a(x) 

>  0 

and  . 

d'(x)  <  0 

then  a 

shock  develops  in  finite  time. 
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Proof :  We  will  prove  part  (3)  only.  The  proofs  for  (1)  and  (2)  are  similar 
v{t,^)  =  x^(t,4).  As  before  characteristics  are  defined  by 

x(o.o.e. 

^  a  a(x)u  u(0,e)  =  , 

when  u(t,e)  =  u(x(  t,4  ) ,  t ).  Let  w(t,e)  = -j|- u(t,C).  Differentiating  (3.2) 
to  5  yields 

|^  =  ^-(u)w  v(0,5)  =  5 

dw 

—  »  a*(x)uv  +  a(x)w  w{0,C>  *  0  • 

Let 


t 

G(t,5)  *  exp  -  /  a(x(s,^))ds  . 

0 


Then  it  follows  from  (3# 3b)  that 


I  w 

w(t,C)  =  g"  (t,5)  /  G(s,€)a'(x(s,C))u(s,5)v(s,5)ds  . 

0 


From  equation  (3.2b)  we  have 


which  implies  that 

u(t,C)  =  "qG”’ (t,C)  . 

(3.6)  combined  with  (3.5)  yields  (the  dependence  on  K  will  be  omitted) 


a 

w(t)  =  G  (t)UQ  /  a'(x(s))v(s)ds  . 
0 


Thus,  the  analysis  of  showing  the  existence  of  a  finite  time  t  at  which  v 
leads  to  the  study  of  the  expression 


Let 


(3.2) 


with  respect 


(3.3) 


(3.4) 


(3.5) 


(3.6) 


(3.7) 


becomes  zero 
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dy 

dt 


^"(u)G“\t)UQ 


t 


I 

0 


a ' (x(s ) ) v(s )ds  . 


(  3.8) 


Let 


so  that  V{0)  =  0  and 


t 

V(t)  =  j  a* {x(s) )v(s)ds 
0 


(3.9) 


V(t)  =  a'(x(t))v(t)  ,  (3.10) 

and  V(  0)  =  a'(C)  from  (3.2).  Hence  the  problem  oC  flndimj  the  time  such  that 

v{T)  =  0  becomes  equivalent  to  showing  V(T)  =  0.  Differentiating  (3.10)  again  we  obtain 

**  ♦  • 

V  =  a'v  +  av 

which  combined  with  (3.9)  and  (3.10)  yields 

a'  •  .  -1 

V  =  —  V  +  aV"G  UpV  , 

V(0)  =  0  ,  (3.11) 

V(0)  =  a’(C)  . 

We  note  that  G  '(t)  >  1  since  a(x)  >  0.  Without  loss  of  generality,  we  can  assume 
that  V(t)  <  0  for  t  S  t0,T],  otherwise  there  exists  T*  >  0  such  that  V(T*)  =  0  and 
by  Rolle's  theorem  there  exists  a  time  T**  such  that  V(T*»>  =0.  With  the  same  argument 
as  above,  we  can  assume  that  V(t)  <  ~S  for  t  >  T, ,  for  some  6  >  0  and  T,  >  u.  Thus 
we  obtain  the  inequality 

V  -  ll  0  -  aV-u^,V  >  0  .  (3.12) 

Multiplying  by  the  integrating  factor  — -  and  integrating  from  0  to  t,  we  obtain 


V - a'(C)  <  u  /  ^"(u(s)  )v(s)ds  , 

a  a  V  s  /  “ 


{  3.  1  3) 


V  *  a'(x(t))[l  +  d  /  (u(s  )  >  V(s  )ds )  .  (1.14) 

0 

Finally,  using  the  hypothesis  on  •fi"  and  the  bound  on  V(t)  we  arrive  at  the  tillowii\g 
estimate  for  V, 

V(t)  >  a'(x(t))(l  -  au^ktl  .  (3.13) 
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Wp  pot-p  tiPrni  in  lirackpts  hpfonpp  nPaativp  in  finite  i  i  ne  p'-rri-ip^  w;*^:  ■'‘• 

sinn  of  a'(x)  ^hows  the  exirtten^e  of  finite  tii»e  at  which  vft)  hecores  .’ere.  T'-i'i 
conpletes  the  proof  of  the  theorem. 

A  similar  resnlt  ran  he  proverl  For  the  erpjation 

+  V  (u)  =  aixK'  •  (n) 

t  X 

U(x^n)  =  ,  (3.  1A) 

where  n^,  as  before,  is  assumed  to  he  a  positive  constant.  On  the  one  hand, the  nonlinear 
dependence  of  a(x)'r'fu)  on  u  causes  the  previous  proof  to  become  more  cumbersome.  Or. 
the  other  hand,  since  the  slope  of  the  characteristics  and  the  nonhomoqeneous  term  have  the 
same  dependence  on  u,  it  is  possible  to  find  u  explicitly  in  terms  of  x  which 
simplifies  the  followinq  proof  considerably. 

Proposition  3.1  Let  ■rMu)  >  0  and  ^^"(u)  >  0. 


r) 

If 

a(x)  <  0, 

,  then  the  solution 

to 

(3.16)  is 

<7lobally  smooth. 

2) 

If 

a(x)  >  0 

and  a' (x)  >  0, 

then 

the  solution  to  (3.16)  is 

qlobally  smooth. 

3) 

Assume  that 

a(x)  A  0,  a’lx)  < 

0 

and  a(x) 

approaches  zero 

as  X  approaches 

infinity.  Suppose  that  there  exists  a  point  5  such  that  /  a(s)ds  <■  ",  Then  a 

a 


Proof. 


shock  develops  in  finite  time. 

Aqain  we  confine  ourselves  to  the  proof  of  part  3.  Refine  characteristics  by 

(u),  x(0,C  )  =  C 


a(x)f>'(u),  u(0,5)  =  Up 


It  follows  from  (3.17)  that 


du 

—  s  a(x) 
dx 


(3.17) 


x(t,C) 

u(t,5)  =  Up  +  /  a(s)dE  . 

K 

Rifferentiatinq  (3.1R)  with  respect  to  5  vield.s 

u^  =  a{x(t,^))x^  -  a(S)  . 


(3.  1S1 


(3.10) 
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On  T-.ne  other  hand/ 


Uj.  can  he  calculated  trom  (3.17^),  i.e.. 


^  dx  _  ^ 

3^  .ft  ■  7t  ’’c 


v-"(u).J^ 


which  combined  With  (3.19)  'jives 


— ! — r  x,^=a(x(t,f))x,-a(f).  'i.2'.  I 

V  (u)  ^ 

Now,  let  f,  tie  the  point  in  the  hypothesis  of  part  3.  Since  ."'(u)  >  )  it  follows  -ta' 
x(t,5)  >  5  for  t  >  0.  (3.18),  then,  implies  that 

|u(t,^)  I  <  Uq  +  M  '  i.  2'  I 


where  M  =  /  a(s)ds.  Thus, 

(u( t,C) )  >  M,  3.22) 

for  some  Since  a(x)  is  positive  and  x(t,C)  >  C  for  t  >  0,  it  follows  that 

u(t,^)  >  Ug.  On  account  ot  <’'Cu)  beinq  positive,  this  in  turn  implies  that 
■;’'(u(t,0)  >  ^'(Uq).  By  (3.17a) 


x(t,C)  >  v-'’(Ug)t  *  C  .  (3.23) 

Since  by  (3.23)  x(t,C)  approaches  infinity  with  t,  there  exists  a  time  T  such  that 
a(x(t,C))  <  ^a(C)  for  t  >  T.  Next  we  show  that  x^(t,C)  <  l  for  c  >  0.  This  follows 
immediately  from  (3.17)  once  we  show  that  u^(t,5)  <  0.  Differentiating  (3.17)  with 
respect  to  5  yields 


du^ 

dt 


a'(x)x^v'(u)  +a(x)o"(u)u^ 


(3.24) 


Since  u(0,C) 


du 

Ug  then  u^(0,C)  =  0.  Theretore  (0,5)  =  a’(C)»"(Ug)  which  is 


negative  by  hypothesis,  i.e.,  u^{t,5)  becomes  negative  tor  a  small  period  ot  time 


du 

the  other  hand,  (3.24)  shows  that  once  u-  is  negative,  — -  remains  negative  due  to  t 

s  dt 

signs  of  a',  •4',  and  The  two  facts  a(x(t,5))  <  ^  att,)  tor  t  ■>  T  and 

x^(t,5)  <  1  enable  us  to  deduce  from  (3.20)  that 

fort>T  (3.2i) 

which  reduces  to 


x^t  *  "  T  t  )  T  . 

Integrating  (3.26)  with  respect  to  t  gives  us 


i  3.  >■ 


-9- 


(■^.27) 


(t)  <  (T)  -  j  -  T)  . 

In<*t|'i>i '  i  ♦  V  (1,?7)  shows  that  t^r  larne  ^nouQh  t  x^(t)  will  herome  netjativp.  We  t^hen 
1.1  to  rompleto  the  proof  ot  this  proposition. 

Sext,  we  tvirn  to  the  ouestmn  of  the  Rieoann  prohlew  for  the  related  eoiiation 

=  n(x,u)  .  (3. 29) 

We  assume  that  V  is  qen'iinely  nonlinear,  i.e.,  >  **  >0.  ConsiHer  (3. 29)  with  the  Piemann 

initial  -onfii  tion 


f  » 

u(x,0)  * 

I 


X  >  0 

X  <  0  . 


(3.2'») 


w.  will  qive  a  brief  outline  of  how  the  local  solution  to  { 3. 2R) -( 3. 2'!)  is  ronstructefl. 
Our  claim  is  *-hat  the  initial  rliscontlnuitv  (3.29)  is  immefiiately  resolver)  by  the 
rorresprinriinrj  conservation  law 

“t  *  ^‘“’x  °  °  •  (3.30) 

Then  the  term  q(x,u)  governs  the  evolution  of  the  resolved  waves.  Hence,  to  solve 
( 3. 2S)-( 3. 29)  we  divide  the  problem  into  two  cases: 

Case  A:  The  solution  to  ( 3, 29)-( 3. 30)  is  a  rarefaction.  Let 

u"(?  )  s 

V  -r  --  - 

be  this  soliitirin,  where 


.  .  . 


b(C  )  if  <  C  <  C 

£  I 

1,  i  f  E  >  f 


(3.31) 


tjj  ”  ”  ■rMu^),  v>'(b(x))  =  X,  ^  ‘ 


^ ons  I 


^  (u),  x(C ,0)  =  0 

o  t 

^  =  ri(x,u),  u(5,0)  =  u*^(5),  C,  <  C  <  5  . 

d  t  A.  r 


(  3.32) 
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It-,  is  not  rhtfioiilr  tn 


Irf*t  (M(r  ,  t )  ,x  (f.  ,  t )  )  ►hp  <5olution  on  ^  ^ 

show  that  '  for  ^  K  ^  K  •  Thus 

.  i  r 

ti(  X,  t )  >=  u(f,  (x,  t ) ,  t) 

Is  a  solution  of  (  3.  30) -( 1 ,  3 1 )  within  the  rpoion  J  <  x  <  x^(t),  with 

X  .  ( t )  =  X  (f.  ,  t ) ,  i  =  r  ,1  • 

1  1 

Case  B;  Ttip  solution  to  ( 3.  2*^ ) -(  3*  30)  is  a  «hoo)<.  Let 


u  (^  )  * 


if  K  >  s 


if  C  >  s 


be  that  solution  with 


•^(u  )  -  «^(u-  ) 
r  I 


Then,  in  a  similar  manner  to  Case  A  we  construct  the  solution  to  ( 3. 28) -( 3. 29) , 


u(x,t)  = 


Uj(x,t)  if  X  <  x(t) 


u^tx,t)  if  x  >  x(t) 


where 


u,  +  '^(u,  )  =  a(x,u.  ) 

1  IX  i 


u . (x, 0)  =  u . 

1  I 


i  =  r,f  , 


and 


~  '^(u^(X,t))  -  '^(Ujj  (x,t) ) 


dx 

at 


x(n)  =  0 


u^(x,t)  -  Uj (x,t) 


(3.33) 


namely. 


(3.34) 


(3.35) 


(3.36) 
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4.  Sohition  try  the  Rif»nann  problem  for  ♦‘he  enuation  (?.!)•  In  '  5 

iterative  ’^rheme  in  onJer  to  obtain  the  solution  to  the  initial  value  :  •'ot  r  *  • 

ueneral  conservation  law 


n  n 

where  y  :  P  ♦  R 
iterative  scheme  is 


+  vMu)^  =  0 
u(x,  0)  =*  '*^(x  ) 

is  smooth  anH  ryenuinely  nonlinear.  T^e  builrtinq 
the  solution  to  the  associated  Riemann  problems 
[  Uj,  X  <  c 


(4.1) 


bloc)i  of  tills 


11  (x,0) 


U  ,  X  >  c  . 

r 


(4.2) 


■rtie  set  of  step  functions  in  (4.2)  is  cliosen  as  a  pointwise  approximation  of  the  initial 
data.  As  shown  in  Section  2,  the  oriqinal  problem  of  the  flow  of  oas  in  a  duct  with  vary- 
ina  cross  section  is  a  two-dimensional  problem.  It  is  the  reduction  of  the  latter  system 
of  eouations  to  a  one-dimensional  Initial  value  problem  that  produces  the  nonhomoqeneous 
terms  in  (2.3).  In  this  section  we  discuss  the  solution  to  the  Riemann  problem  for  (2.3) 
which  arises  from  the  discretization  of  the  Initial  data  and  the  boundary  of  the  duct.  It 
is  this  simultaneous  discretization  that  maltes  our  treatment  of  the  Riemann  problem  differ¬ 
ent  from  the  one  discussed  in  (51.  First  we  note  that  (2.3)  can  be  written  in  the  form 


(Ap  )^  +  (Apu)  =  0  , 

t  X 

“t  *1  '“"’x 


(4.3) 


where  P  is  defined  by 


P(P) 


P 

/ 


P’  (s) 
s 


ds 


and  A(x)  has  the  form 


(4.4) 


A(x) 


1  , 

1  -  e. 


X  <  O 
X  >  0 


(4.5) 
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I 

i  (P  ,u  ) 

^  4-  + 


y  <  0 
X  >  0  • 


(4.6) 


As  in  f5),  we  assume  that  both  the  initial  con^iition  and  the  boundary  of  the  duct 
have  small  hounded  variations.  When  e  =  o  (4.3)-<4.6)  reduces  to  the  classical  Riemann 
problem  for  the  eauations  of  qas  dynamics  in  a  uniform  tube  (81.  For  the  case  e  positive 
we  apply  the  same  ideas  as  in  Section  3,  namely,  the  solution  to  (4.3)-(4.6)  can  be  viewed 
as  a  small  perturbation  of  the  solution  to  the  correspondino  problem  when  e  *=  o.  The 
implicit  function  theorem  is  the  main  tool  in  obtaininq  the  exact  solution  of  {4.3)- 

(4.6).  To  illustrate  the  method  we  choose  a  particular  solution  of  the  e  »  0  case  and 
carry  out  the  necessary  calculations.  Let 


(P_,u_), 

0  <  -^  < 

s 

(P_.u  ), 

t 

X 

s  <  -  < 

(p’(p 

m 

q(f),  (P’(P  (p>{p^)/''2 

t  m  t  + 


be  the  physically  admissible  solution  to  (4.3)-(4.6)  with  e  =  Or  i.e.,  the  solution  to 

the  Riemann  problem  consists  of  a  backward  shock  (p  ,u  ;  p  ,u  ;  s)  and  a  forward 

-  -  m  m 

rarefaction  wave  connectinq  (p  ,u  )  to  (p  rU  )  (cf.  (81).  Then  che  solution  to  (4.3)- 

mm  +  + 

(4.6)  with  e  positive  consists  of  a  backward  shock  (p  ,u  ;  p^(e),u^(e);  s(b)),  a 
discontinuity  (P  ^  (e  ),u^  (e  ) ;  0  )  ,n^(.e  ) ;  O)  which  is  due  to  the  oeometry  of  the  duct,  and 

a  forward  rarefaction  connectinc  (p  ^(e  ),vi  ^(e  ) )  to  (p^,u^).  The  five  formulae  relating 
s(E),  p^(E),  u^(E),  p^(E),  and  U2(e  )  are 


s(p , 


p_)  =  P,u^  -  P_u_ 


2  2 

s  (P  ^u^  -  P  Jij  =  P  ^>1^'  +  P(P  , )  -  P_u_  -  p(P  _)  , 
(1  *  ^  ^P2“2  ‘’l“l  ’ 


(4.7) 


u,  *  / 

P 


p'  (P  ) 


do  , 
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In  that  r:ase  we 


Equations  (4.7)  are  simplified  considerably  after  eliminatinn  s  and  u^. 
can  formulate  the  above  problem  in  the  form 


where  ^  • ''i '  ^2 ' ^3 ' 


(4.8) 


P, 

F,  (P,  ,P  2<“2'®  *  °  ~  ^  2^2  ■  *’  l“-  *  1*  -  p(P_)  )  (P  ,  -  p_)]'^  , 


2'%'‘’2'“2'‘^'  “  "2  -  u+  -/ 


/p'  (P  ) 


dp  , 


(4.9) 


We  point  out  that  the  aiqn  of  the  square  root  in  (4.9)  is  chosen  so  that  the  usual  entropy 
condition  is  satisfied  (1).  (4.8)-(4.9)  is  now  set  up  for  applyinq  the  implicit  function 

theorem,  the  problem  is  solved  if  we  can  uniquely  determine  p^,  p^.  and  Uj  in  terns  of 
e,  Tto  this  end  we  calculate  the  Jacobian  of  F  with  respect  to  P^.  P^i  and  Uj  at 
e  •  0.  It  can  be  shown  that  (all  partial  derivatives  are  evaluated  at  P^  »  P  ,  P^  *  P  , 
Uj  «  u_,  and  e  =  0) 


1  rP  . 


=  -U  +  (p'  (P  )  ) 


1/ 


■1#P. 


1  »u. 


2,p, 


0, 


F 


2rP 


2 


(p'  (P  ))  ’*^2 


2rU, 


^■I  a  =  P  U^  -  P  p'  (P  )  .  P,  -  =  -P  U^  +  P  P’  (P  ).  F  =  0  . 

3,p,  -  -  -  -  3,P^  -  -  -  -  S.u^ 

A  simple  calculation  shows  that 


det 


3(Fi,F2,P3) 


=  2P  (u^  -  p'(P  ))(P’(P  ))^^  . 


E«0 


(4.10) 


Thus,  if  the  initial  step  (P  ,u  ;  P^,u^)  is  such  that 

-  P'(P_)  (4.11  ) 

i.e.,  the  Upstream  flow  of  the  qas  is  either  subsonic  or  supersonic,  we  see  that  the 
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implicit  function  theorem  can  he  applieo  to  solve  uniquely  for  P^,  P and  U2  in  terms 
of  e.  This  completes  the  solution  to  the  Riemann  problem  (4. 

There  are  still  two  interestinq  problems  in  connection  with  (4,3)  whose  answers  would 
be  valuable  both  to  the  theory  and  the  application.  The  first  question  is  whether  the 
above  scheme  actually  converqes  to  the  weak  solution  of  the  initial  value  problem.  The 
second  question  in  how  easily  this  scheme  can  he  implemented  numerically.  We  are  presently 
studyinq  these  questions. 
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